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Abstract 

Rankine-Hugoniot condition has been solved to study phase transition in 
astrophysical scenario mainly in the case of phase transition from neutron 
star (NS) to quark star (QS). The phase transition is brought about by a 
combustion front, which travels from the center to the surface. The equations 
of state and temperature plays a huge role in determining the nature of the 
front propagation, which brings about the phase transition in neutron stars 
(NSs). Magnetic field has been introduced and the modified conservation 
condition for the perpendicular and oblique shocks is obtained. Numerical 
solution of the perpendicular shock has been shown in the figures, which 
finds that the magnetic field helps in shock generation. It indirectly hints at 
the instability of the matter and thereby the NS for very high magnetic field, 
implying that NSs can only support finite magnetic field strength. 
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1. Introduction 

When the velocity of a fluid in motion becomes comparable with or ex- 
ceeds that of the sound, the effect due to compressibihty of the fluid become 
of prime importance. For a wave propagating in a non conducting gas, when 
the amplitude is so small that linear theory applies, the disturbance prop- 
agates as a sound wave. If the gas has a uniform pressure and density, the 
speed of propagation of sound and the wave profile maintains a fixed shape, 
since each part of the wave moves with same speed. However, when the wave 
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possesses a finite amplitude, so that nonlinear terms in the equation become 
important, the crest of the sound wave moves faster than its leading or trail- 
ing edge. This causes a progressive steepening of the front portion of the 
wave as the crest catches up and ultimately, the gradient of pressure, den- 
sity, temperature and velocity become large that dissipative processes, such 
as viscosity or thermal conduction are no longer negligible. Then a steady 
wave shape is attained, called a shock wave. The shock wave moves at a 
speed in excess of the sound speed, so the information cannot be propagated 
ahead to signal its imminent arrival, since such information would travel at 
sound speed relative to the undisturbed medium ahead of the shock. The 
dissipation inside the shock front leads to a gradual conversion of the energy 
being carried by the wave into heat. Thus, the effect of the passage of a 
shock wave are to convert ordered (flow) energy into random (thermal) en- 
ergy through particle collisions and also to compress and heat the gas. The 
shock front itself is in reality a very thin transition region. Its width is typi- 
cally only a few mean-free paths, with particle collisions establishing the new 
uniform state behind the shock. 

The relativistic shock propagates into a medium with a changing equation 
of state. Therefore, a simple analysis of the jump condition for a polytropic 
or perfect fluid is not adequate and a deep understanding of this problem calls 
on for the full theoretical description of the relativistic shock in a medium 
with arbitrary equation of state. Further complication might arise if there 
is a presence of significant magnetic field. In fact, one can show that the 
relative importance of a magnetic field can grow during a collapse. In the 
field of nuclear physics, high energy collisions among heavy ions can be mod- 
eled by using fluid dynamical concepts. Also, some current models under 
investigation predict that relativistic shocks (or relativistic detonation and 
deflagration) might be related to the phase transition from nuclear matter 
to a quark matter. 

Relativistic shock waves have been the subject of early investigation in 
relativistic fluid dynamics and magneto-fluid dynamics. In relativistic fluid 
dynamics the pioneering work is that of Taub [l| where the relativistic form 
of jump condition is established. A detailed analysis of the thermodynamic 
treatment of classical shock waves, is due to Thorne ^ . Explicit solutions of 
the jump condition have ben obtained for special equation of states. Shock 
wave in relativistic magneto-fluid have been investigated extensively and in 
rigorous mathematical way by Lichnerowicz. Detonation and deflagra- 
tion waves in relativistic magneto-fluid dynamics for nuclear physics and 
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cosmology have been investigated by 0, [sj . 

In Astrophysical scenario shock plays very important role in determining 
the outcome of compact stars. In the case of massive stars, in the range 
between 8-100 solar masses, which are thought to be the progenitors of type 
II supernovas, one of the most viable mechanism for producing a explosion is 
gravitational collapse and bounce In this case a shock is formed outside 
the inner core, which propagates outwards reaching relativistic speeds. Shock 
waves are also responsible for phase transition and gamma ray bursts (GRB) 
in compact stars. 

In this paper I will mostly concentrate on the phase transition scenario 
in a compact star. Witten [7j] conjectured of strange quark matter (SQM), 
consisting of approximately equal numbers of up (u), down (d) and strange 
(s) quarks, is believed to be the ground state of strong interaction. This was 
supported by model calculations for certain ranges of values for strange quark 
mass and strong coupling constant After that there has been constant 
efforts at confirming the existence of and SQM, though transiently, in ultra 
relativistic collisions. On the other hand, SQM could naturally occur in the 
cores of compact stars, where central densities are expected to be an order of 
magnitude higher than the nuclear matter saturation density. Thus, neutron 
stars which have sufficient high central densities might convert to strange 
star, or at least hybrid (a star with a quark core) stars. These transitions may 
lead to various observable signatures in the form of a jump in the breaking 
index and gamma ray bursts jo], [lo[ , and a full QS might help in explaining 
the phenomena of observed quasi periodic oscillations [1 Ij . 

There may be several scenarios by which neutron stars could convert 
to quark stars. It may happen through a "seed" of external SQM 12|], or 



triggered by the rise in the central density due to a sudden spin-down in older 



neutron stars 13 . Several authors have studied the conversion of nuclear 



matter to strange matter under different assumptions [ij, ll5|, ll6|, ll7|, [18 



19, 
25| 



20l . l21l |22| . |23| . l24j | . They have been summarized in a recent work of ours 
and for the constraint of space, I do not repeat them here. 

After the discovery magnetars, some compact stars were found to have 
very high surface magnetic fields. So the study of relativistic rankine-hugoniot 
condition is not sufficient. To have a full understanding of the properties of 
NS and its phase transition to QS, such conditions should be examined in 
the presence of high magnetic fields. In this paper I wish to carry out such 
a basic calculation keeping our focus mainly on the astrophysical scenario. 
The paper is organized as followed: first I will discuss the general rankine 
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hugoniot condition as a discontinuity in the conversion front. In section III, 
I will introduce magnetic field and study the new set of modified conserva- 
tion equations. Then in section IV, I shall show my results and finally I will 
discuss and summarize them. 



2. General Rankine-Hugoniot condition 



In this section I will discuss the general rankine-hugoniot condition as the 
conservation equations which balances the conversion of neutron proton (n-p) 
matter to two-flavour quark matter, consisting of u and d quarks along with 
electrons for ensuring charge neutrality. I heuristically assume the existence 
of a combustive phase transition front. Using the macroscopic conservation 
conditions, I examine the range of densities for which such a combustion 
front exists. 

Let us consider the physical situation where a combustion front has been 
generated in the core of the Neutron star. This front propagates outwards 
through the neutron star with a certain velocity, leaving behind a u-d-e mat- 
ter. In the following, I denote all the physical quantities in the hadronic sector 
by subscript 1 and those in the quark sector by subscript 2. The conservation 
condition for energy-momentum and baryon number relates the quantities on 
the opposite sides of the front. In the rest frame of the combustion front, 
these conservation conditions is given by 2J, |26|, l27l |: 



^i^i7i + Pi 



^2^272 +P2, 



(1) 



and 



(2) 



(3) 



In the above three conditions f , (i= 
— ^ is the Lorentz factor, oji = 



=1, 2) is the velocity, Pi is the pressure, 
■ ei + Pi is the specific enthalpy and 



is the energy density of the respective phases. 

The velocities of the matter in the two phases, given by equations ([I]|3]), 



can be solved, such that 26 



(P2 -Pi)(e2 

(£2 - ei)(ei + P2) ' 



(4) 
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and 

2 _ {P2 - Pi){ei + P2) 

It is possible to classify the various conversion mechanisms by comparing 
the velocities of the respective phases with the corresponding velocities of 
sound, denoted by Csi, in these phases. Thes conditions are summarized in 

3- 

For the conversion to be physically possible, velocities should satisfy an 
additional condition, namely, < < 1. Here I find that the velocity 
condition puts severe constraint on the allowed equations of state. 



3. Magnetic field Inclusion 

In a conducting g gnetic field can interact strongly with the flow. 

The analysis of the shock waves therefore becomes more complex, but the 
basic principles remains the same. A set of jump condition can again be 
derived, but they are considerably more complicated than the pure hydro- 
dynamic shock case. The extra complexity arises both from the presence of 
extra variable, namely the magnetic field strength, and also from the fact 
that the magnetic field and the matter velocities may be inclined with the 
shock normal. 

A shock propagating through an magnetic fluid produces a significant dif- 
ference in matter properties on either side of the shock front. The thickness 
of the front is determined by a balance between convective and dissipative 
effects. However, dissipative effects at high temperature are only comparable 
to convective effects when the spatial gradients in matter variables become 
extremely large. Hence, shocks in such matter tend to be extremely narrow, 
and are well-approximated as discontinuity. The hydrodynamical equations, 
and Maxwell's equations, can be integrated across a shock to give a set of 
jump conditions which relate matter properties on each side of the shock 
front. If the shock is sufficiently narrow then these relations become inde- 
pendent of its detailed structure. 

In the rest frame of the shock, the shock front coincides with the y-z 
plane. Furthermore, the regions of the plasma upstream and downstream 
of the shock, which are termed regions 1 and 2, respectively, be spatially 
uniform and non-time- varying. It follows that d/dt = d/dy = d/dz = 0. 
Moreover, d/dx = 0, except in the immediate vicinity of the shock. Finally, 
let the velocity and magnetic fields upstream and downstream of the shock 
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all lie in the x-y plane. The magnetic field is given by Bi for the respective 
phases. 

The first nontrivial shock is called perpendicular shock in which both the 
upstream and downstream plasma fiows are perpendicular to the magnetic 
field, as well as the shock front. The conservation condition are given by 

B ^ B 
uJivf-fl +Pi + -^ = ^2vhi + ^2 + , (6) 

B ^ B ^ 

uJiVijl + ^171-— = a;2V272 + ^^272-;—, (7) 
47r 47r 

BlVi'Ji = 52^272 (8) 

and 

niVi^i = n2V2^2- (9) 
The first three equation can be reduced to two equation, given by 

2 2 I I -^1^ 2 2 I I -Si^/'^^171n2 /■,n\ 

OTT OTT f272 

UiVi'j^ + ^^171-;— = ^^2^^272 + ''^272-; , (H) 

47r 47r V272 

I now solve for Vi nad V2. 

The most general shock is the oblique shock in which the plasma velocities 
and the magnetic fields on each side of the shock are neither parallel nor 
perpendicular to the shock front. The Rankine-Hugoniot condition are given 

by 



2 2 I I -^1 -^la: 2 2, , -°2 ^2x 
^iV^x^l +Pl + ^ = ^2^2^72 + P2 + ^ , (12) 



i^lVlxVlyJl + Vi^-fi = UJ2V2xV2y^2 + ^2x72 Z , (13) 

47r 47r 



2 , Bi' B,,(Bi.Vi)^i 2 , B2' B2x{B2.V2h'2 
^^iVlxll + Vlxll- = l^2V2xl2 + V2xl2-, , 

47r 47r 47r 47r 

(14) 
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BiyVixIl - VlyBi^-fi = B2yV2x'l2 - V2yB2xl2 (15) 



Bix = B2x (16) 

and 

niVix'Ji = n2V2x72- (17) 

There may be two easier cases for the above comphcate equation. 
Case 1. viy = = V2y 

Then the first four equation would simphfy to 

B ^ B ^ 

^i^Lt? +Pi + -^= ^2vl^-il + P2 + (18) 

2 B^^ ^i.(^i.t;i)7i 2^ ^2 BUB2-V2h2 .... 

U:,Va,+V,l,— =UJ2V2l2+V2l2— , (19) 

B2y = ^B,y (20) 

^272 

Case 2. B^y = = B2y 

^l^lxll +Pl= t^2f 2x72 + P2, (21) 



2 , Bi,(Ei.Wi)7i 2 B2' BixiB2.V2h2 
^^1^^1x71 + ^1x71-^ = ^2V2xl2 + ^2x12- , 

(22) 

^^iVlxVlyJl'^ = W2^^2xf2y72^ (23) 

4. Results 

I start my calculations by using the nuclear matter EOS obtained through 



nonlinear Walecka model [29[. In the present paper, I consider the conver- 
sion of nuclear matter, consisting of only nucleons, to a two-flavour quark 
matter. The final composition of the quark matter is determined from the 
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nuclear matter EOS by enforcing the baryon number conservation during the 
conversion process. While describing the state of matter for the quark phase 
1 consider a range of values for the bag constant. Nuclear matter EOS is 
calculated at zero temperature, whereas, the two-flavour quark matter EOS 
is obtained both at zero temperature as well as at finite temperatures as 
during the propagation of the shock it may heat up the matter. 

To examine the nature of the hydrodynamical front, arising from the 
neutron to two-flavour quark matter conversion, 1 plot, in flg.l, the quanti- 
ties vi,V2,Csi and Cs2 as functions of the baryon number density (ns). As 
mentioned earlier, the u and d quark content in the quark phase is kept 
same as the one corresponding to the quark content of the nucleons in the 
hadronic phase. With these fixed densities of the massless u and d quarks 
and electrons, the EOS of the two-flavour matter has been evaluated using 
the bag model prescription. 1 find that the velocity condition ( > ) is 
satisfied only for a small window of ±5.0MeV around the bag pressure 
S^/^ = 160 MeV. The constraint imposed by the above conditions results in 
the possibility of defiagration, detonation or supersonic front as shown in the 
fig. 1. 

In fig. 1, 1 considered both the phases to be at zero temperature. A 
possibility, however, exist that a part of the internal energy is converted 
to heat energy, thereby increasing the temperature of the two-flavour quark 
matter during the exothermic combustive conversion process. In fig. 2, I plot 
the variation of velocities with density at T = 50MeV, for which significant 
change is noticed. This figure show that the range of values of baryon density, 
for which the flow velocities are physical, increases with temperature. In the 
present paper I have considered only the zero temperature nuclear matter 
EOS. On the other hand, equation of state of quark matter has a finite 
temperature dependence and hence the difference between vi and V2, varies 
with temperature. 

The variation of the velocities with temperature is due to the fact that 
higher temperature means higher energy. As the energy increases the parti- 
cles becomes more energetic which means the matter becomes more excited 
and compressible. As matter becomes more compressible due to increase in 
temperature, now there is a chance for shock formation which was previously 
not possible. 

Till now I have used equation for the normal hydrodynamics for the gen- 
eration of the curves. The change in temperature in the EOS accounted for 
the change in the nature of the graph. Now I will plot graphs solving mag- 
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Figure 1: Variation of different velocities with baryon number density for T = MeV, 
B^/'^ = l&QMeV and strange quark mass = 2QQMeV. The dark-shaded region cor- 
respond to deflagration, light-shaded region correspond to detonation and the unshaded 
region correspond to supersonic conversion processes. 




number density 



Figure 2: Variation of velocities with baryon number density for T = 80 MeV B^/'^ = 
IQOMeV and strange quark mass = 200MeV. Different regions correspond to different 
modes of conversion similar to that of fig 1 . 
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netic field induced hydrodynamics. I solve the new conservation conditions 
for the perpendicular shock in which both the upstream and downstream 
plasma flows are perpendicular to the magnetic field, as well as the shock 
front. 

In fig. 3 I have plotted curves for different velocities with baryon density 
for the zero temperature case. The magnetic field used for the generation 
of the curve is i? = 5 x lO^^G. I find that due to introduction of the mag- 
netic field the range of values of baryon density, for which the flow velocities 
are physical, increases. Its nature is quite similar to that for the case of 
temperature. Due to the introduction of the magnetic field the the pressure 
increases (pressure due to magnetic field is B'^/Sn) for the same value of 
baryon number density. The magnetic field in the respective phases adjust 
in such a way that the resultant energy and pressure of the both phases gets 
modified to make the range of baryon number density to increase. For both 
higher and lower values of baryon density, previously there was less chance of 
shock formation, but now due to the introduction of the magnetic field there 
is a greater chance of shock formation. And therefore the range of baryon 
density gets much wider. 

Next I plot curves for different velocities for the finite temperature case 
T — SOMeV. For same value of the magnetic field B — 5x lO^^G, the range 
of baryon density gets much more wider. This is due to the fact that now 
both the temperature and the magnetic field work hand by hand to increase 
the chances of shock generation. Both process ensures by its own way that 
the matter parameters adjust itself in such a way that there is a greater 
probability of shock generation. 

In fig. 5 1 have plotted exclusively vi and V2 for two cases, one without 
magnetic field and the other with magnetic field. The nature of the curve 
remains same, that is Vi is always greater than V2, which means the shock 
front propagates outward of the star. The range of baryon density, for which 
the matter velocities are finite, increases with magnetic field. I have plotted 
this for the zero temperature case. In fig. 6 I have plotted the same for the 
finite temperature case and find that for lesser value of magnetic field, same 
increase in range of baryon number density is seen. All the reasons for this 
nature is explained in the previous paragraph. 

In fig 7. I have plotted vi for different values of magnetic field. I find 
that as the value of the magnetic field increases the range of baryon number 
density increases which is what was expected (and discussed previously). 
But as we go on increasing the magnetic field Vi does not comes down on 
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Figure 3: Variation of different flow velocities with baryon number density for T = MeV 
and magnetic field of S = 5 x lO^^G. 

the lower half of the curve. This is due to the fact that, at such high value 
of the magnetic field the matter becomes unstable. At other half where 
matter is at much higher density, it can support such field strength. But, if 
we further increase the magnetic field the matter cannot support such high 
fields whatever the density might be. The maximum value of magnetic field 
that matter can support is few times lO^^G (in our case the cut off value is 
2 X 10^ ''(?). So I find that there is a cut off value for the magnetic field, and 
from here we can indirectly say that NS also can support up to a finite value 
of magnetic field. This curve is plotted for zero temperature, and in fig. 8 
I have plotted the same for finite temperature. Qualitatively the nature of 
the graph remains the same only the quantitative value of the magnetic field 
changes. It shows that hotter matter can support lesser value of magnetic 
field than the colder one. 

5. Summary and discussion 

Finally in this section 1 summarize my results. 1 find that rankine- 
hugoniot condition can be solved to determine the condition for different 
types of wave generation in a neutron star. It also determines the mode of 
the propagation of the wave front. The temperature (finite temperature of 
the matter) helps in the generation of the front and as the temperature rises 
the wave front can generate both at much lower and at much higher baryon 
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Figure 4: Variation of velocities with baryon number density for T = 9S)MeV, and mag- 
netic field of B = 5 X lO^^G. 
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Figure 5: Comparison of Vn and Vs with baryon number density for T = OMeV for two 
values of magnetic field strength B = OG and B = 1 x lO^^G. 
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Figure 6: Comparison of w„ and Vs with baryon number density for T = SOMeV for two 
values of magnetic field strength B = OG and B = 5 x lO^^G. 
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Figure 7: Comparison of v„ with baryon number density for T = OMeV for different 
values of magnetic field strength. 
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Figure 8: Comparison of with baryon number density for T = %QMeV for different 
values of magnetic field strength. 

densities which was not possible for the zero temperature case. Next I write 
down the modified conservation conditions in presence of magnetic field. I 
have written down conditions for both the perpendicular and oblique shock 
waves. The inclusion of the magnetic field introduces not only extra condi- 
tions but also the earher existing conditions gets modified. I have solved the 
velocity of the matter of the two phases and plotted curves for the perpen- 
dicular wave, which was not obtained analytically. The conditions cannot 
be solved analytically and and therefore 1 have solved the nonlinear simul- 
taneously equation numerically. I have matched my nonmagnetic numerical 
results with the analytically solvable nonmagnetic solutions. The oblique 
wave equation gets very much complicated and the simultaneous equations 
does not converges for different value of the baryon densities. So I have not 
plotted the results of the oblique waves. 

Solving the perpendicular wave for finite magnetic field I find that the 
range of baryon density, for which the fiow velocities of matter are physical, 
increases with magnetic field strength. This is because by the introduction 
of the magnetic field the resultant pressure and energy redistribute in such a 
way that, for same baryon density, the changes for shock generation increases. 
I also find that there is a cut off magnetic field strength that matter can 
support, beyond which the matter becomes unstable and the ffow velocities 
becomes imaginary. This on the other hand suggest that a NS cannot support 
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magnetic field beyond a certain field strength. To finally summarize our result 
I mention that this is the first instance where such a treatment of modified 
conservation condition has been done in the presence of magnetic field in the 
astrophysical phase transition scenario. This provides with new interesting 
results, that has not been anticipated before and also indirectly hints at the 
instabihty of the NS at very high magnetic field. More interesting results is 
anticipated if the full oblique wave equations can be solved and I am now 
focussing mainly on that problem. 

I would hke to thank Grant No. SR/S2HEP12/2007, funded by DST, 
India for financial support. 
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